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Abstract 

We present universal factorization formulas describing the behavior of one-loop QCD 
amplitudes as external momenta become either soft or collinear. Our results are valid to all 
orders in the dimensional rcgularization parameter, e. Terms through C(e 2 ) can contribute 
in infrared divergent phase space integrals associated with next-to-next-to-leading order jet 
cross-sections. 



1 Introduction 



In recent years, significant progress has been made in computing the next-to-leading order 
(NLO) corrections to multi-jet rates within perturbative QCD. This progress has come in the 
form of both matrix element calculations (see e.g. refs. |l|, |j) and algorithms for numerical jet 
calculations (see e.g. refs. H ||). An important next step would be the computation of next- 
to-next-to-leading-order (NNLO) corrections to multi-jet rates. As an example, a calculation of 
NNLO contributions to e + e~ — ► 3 jets is needed to further reduce the theoretical uncertainty in 
the determination of the QCD coupling, a s , from event shape variables ||. In fact, many of the 
experimental searches for new physics require a precise understanding of the QCD background, 
and a more accurate determination of a s . At hadron colliders, NNLO calculations would reduce 
the dependence of the multi-jet rates on the factorization and renormalization scales, and would 
allow a more detailed study of jet structure. 

To compute n particle production at NNLO, three sets of amplitudes are required: a) n 
particle production amplitudes at tree level, one loop and two loops; b) n + l particle production 
amplitudes at tree level and one loop; c) n + 2 particle production amplitudes at tree level. 
For example, the computation of NNLO e + e~ — > 3 jet production requires e + e~ —* 5 parton 
amplitudes at tree level ||, e + e~ — > 4 parton amplitudes at tree level J7| and at one loop |8], 
and e + e~ — > 3 parton amplitudes at tree level [{Tcfl , one loop [[H]] and two loops. In terms of the 
amplitudes required, the crucial missing piece is the two-loop calculation. In fact, no two-loop 
computations exist for configurations involving more than a single kinematic variable, except in 
cases of maximal supersymmetry (T^j. 

In general, loop-level amplitudes have a complicated analytic structure, with the complexity 
increasing rapidly with the multiplicity of kinematic variables. The analytic structure simplifies 
in the infrared limit (where parton momenta are soft or collinear) when the amplitudes become 
singular. Infrared singularities exhibit universal, i.e. process independent, behavior, manifesting 
themselves as poles in the dimensional regulator e = (4 — D)/2 after integration over virtual or 
unresolved momenta. The Kinoshita-Lee-Nauenberg theorem [13| guarantees that the infrared 
singularities must cancel for sufficiently inclusive physical quantities when the real and virtual 
contributions are combined. For processes without strongly interacting initial state particles, 
like e + e~ — > X, the cancellation is complete. For processes with strongly-interacting initial 
states, like jet production in hadron collisions, initial-state collinear divergences survive the 
cancellation and are factorized into the parton distribution functions, reducing the dependence 
of the hadron cross section on the factorization scale, up. 

The structure of infrared singularities at NLO is well understood. In the squared tree- 
level amplitudes, any one of the n + l produced particles can be unresolved in the final state. 
The ensuing infrared singularities are accounted for by tree-level soft (14, ^] and collinear 
splitting functions 16]. These have also been combined into a single function [17]. In addition, 
the universal structure of the coefficients of the 1/e 2 and 1/e poles for the one-loop virtual 
contributions to n-particle productions is known (||, 18, 19]. The universality of the structure 
of the singularities has been exploited in building general-purpose algorithms |, |, 01 for NL0 
jet production calculations. 

The study of the infrared structure at NNLO is still underway. A detailed understanding of 
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the infrared singularities that arise from the combination of virtual loops and unresolved real 
emission will be crucial to the development of methods for performing NNLO calculations. In 
the squared tree-level amplitudes, any two of the n + 2 produced particles can be unresolved. 
The resulting soft, collinear, and mixed soft/collinear singularities are described by universal 
tree-level double-soft [[0|, double-splitting and soft-splitting [2C] functions, respectively. In 
addition, the universal structure of the coefficients of the 1/e 4 , 1/e 3 and 1/e 2 poles has been 
determined [21] for the two-loop virtual contributions to n-particle production. 

In the interference between a one-loop amplitude for n + 1 particle production and its tree- 
level counterpart, unresolved real emission generates additional 1/e 2 divergences, which when 
combined with one-loop virtual singularities brings the total order of divergence to 1/e 4 . In 
order to evaluate the interference terms to O(e ), the (n+ l)-parton one-loop amplitude must be 
evaluated to C(e 2 ). In light of the already complicated analytic structure of one-loop amplitudes 



to O(e ) [||, ||, 22, ^3], 24, 25 1, such a calculation would be a truly formidable task. A more 
reasonable approach is to evaluate the amplitudes to higher order in e only in the infrared 
regions of phase-space where the amplitudes factorize into sums of products of n-parton final- 
state amplitudes multiplied by soft or collinear splitting amplitudes. It is only these soft and 
collinear splitting amplitudes and the n-parton final-state one-loop amplitudes that must be 
evaluated to higher order in e. 

In this paper we describe in detail the methods used to obtain the soft and collinear splitting 
amplitudes to all orders in e and present new results for the soft and splitting amplitudes with 
quarks, which have previously been given only through O(e ) p4| , |26|| . We present our results 
both in terms of the helicity representation and in terms of formal polarizations and spinors. The 



latter form is convenient if one is working entirely in conventional dimensional regularization [ 27 1 . 



In a previous letter [28], we presented all orders in e results for the one- loop pure gluon soft and 



collinear splitting amplitudes and used the one-loop soft amplitudes to 0(e) to re-derive [29] 



next-to- leading logarithmic corrections to the Lipatov vertex [ 30 1 . 

In order to present the factorization properties in their simplest forms, we decompose one- 
loop QCD amplitudes first into partial amplitudes [31, 32|, which follow from the color structure, 



and then into primitive amplitudes [||, E4j. The benefit of this last step is that color factors are 
completely disentangled from primitive amplitudes, allowing for relatively simple description of 
the soft behavior of the amplitudes. 

This paper is organized as follows: In section |2|, we review the decomposition of one-loop 
QCD amplitudes by color ordering and the further decomposition into primitive amplitudes; 
in section |||, we review the collinear and soft factorization properties of QCD amplitudes; in 
section we describe our method for obtaining the soft and collinear splitting amplitudes to all 
orders in e, which is based upon the discussion of ref. [33|, and present our results; in section [5], 



we check our result by comparing to those obtained from N = 4 supersymmetric amplitudes 
(which are known to all orders in e ]p4| ) and those obtained from gggH amplitudes using the 
effective ggH coupling [35] generated by a heavy fermion loop, again being careful to keep all 



higher order in e contributions. (For a discussion of this calculation through C(e°) see ref. p(J.) 
We note that an interesting technique for obtaining the collinear splitting amplitudes has 



been recently proposed |p7| , using the unitarity reconstruction methods of refs. [^, 38 1. This 
method appears to allow for a straightforward generalization to higher loops. 
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2 Review of color decompositions and primitive amplitudes 



To discuss the properties of one-loop QCD amplitudes as momenta become either soft or collinear 
we will use 'primitive amplitudes'. Primitive amplitudes, defined in refs. fi ^4|] and originally 
motivated by the structure of fermionic string theory amplitudes, are gauge invariant building 
blocks from which QCD amplitudes, including their color factors, can be built. An important 
characteristic of primitive amplitudes is that they have a fixed ordering of the external legs, lead- 
ing to a relatively simple analytic structure. In traditional representations of QCD amplitudes, 
factorization in the collinear limit is fairly straightforward, but in the soft limit, color factors 
become entangled with kinematic factors in a non-trivial way. At tree level, a conventional color 



decomposition [ 31 1 is sufficient to disentangle the soft factorization properties. This is no longer 
true at one-loop where the presence of sub-leading color structures results in an even deeper 
entanglement of color and kinematics. Our purpose in using primitive amplitudes is to provide 
a clean factorization of one-loop amplitudes in the soft and collinear regions and to separate 



color issues from the kinematic issues. In section 4.6.3 we will make use of the formulas reviewed 
in this section to give an example of how the simple soft factorization of primitive amplitudes 
involving quarks leads to a nontrivial tangle in terms of the more conventional color ordered 
partial amplitudes. In ref. |28|], where only the pure glue case was dealt with, there was no need 
to introduce primitive amplitudes since the standard leading color partial amplitudes play the 
same role. 

For the cases of ra-gluon amplitudes and two quark, n — 2 gluon amplitudes general for- 
mulae expressing color ordered partial amplitudes in terms of primitive amplitudes have been 



presented |24|, [2(|. For the purposes of this paper we review this one-loop decomposition for 
the limited cases of four or five colored partons, corresponding to NNLO pp — > 2 jets, which 
will likely be among the first cases to which the results in this paper can be applied. We also 
very briefly outline the construction of the primitive decomposition of four quark, n — 4 gluon 
amplitudes. The primitive decomposition for the one-loop contributions to e + e~ — > 3 jets and 
e + e~ — ► 4 jets may be found in ref. ||]. 

The usual color decomposition follows from replacing the color structure constants with 
commutators of fundamental representation matrices,^ 

jabc = — L(TY[T a T fe r c ] - Ti[T b T a T c }) , (2.1) 
v2 

and then using SU(N C ) Fierz identities, 

to combine traces. In the case of amplitudes with purely adjoint particles one may use instead 
the U(N C ) Fierz identities which do not contain the second self contraction term. 

At tree level the pure gluon amplitude may be expressed as (I], |3l| ] 

A* ee (l,2,...,n) = g n - 2 Tr(T a °WT a °w • • ■ TM«)) ^ ree (tr(l), a(2), ... ,a(n)), (2.3) 

n 



*Note that we use a non-standard normalization for fundamental representation matrices, Tx[T a T b ] = S ab 
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where S n /Z n is the set of all non-cyclic permutations. For the case of two quarks in the funda- 
mental representation and (n — 2) gluons in the adjoint representation the color decomposition 
is 

^T(l g -,2 g ,3,...,n) = «f- 2 E (TM3)...r<M^ ^(1,-, 2 q ; <r(3), . . . , a(n)) , (2.4) 



ires'; 



n-2 



where S n -2 is the permutation group on n — 2 elements. Similar expressions exist for the case 
of larger numbers of fermions. In each case, the terms which multiply different color structures 
are called 'partial amplitudes'. 

We may also define a set of tree-level primitive amplitudes for the pure gluon case to corre- 
spond to the partial amplitudes. For the case of two fermions the primitive amplitudes are the 
set of partial amplitudes one would obtain if all particles, including the fermions, were in the 
adjoint representation. 

This decomposition yields color structures consisting of products of fundamental represen- 
tation matrices. In the simplest color structures the product involves a single chain of repre- 
sentation matrices. Tree-level amplitudes only use the leading color structures. There are also 
more complicated (sub-leading) color structures, which come in one-loop, in which the product 
breaks into more than one chain. 

For example, the color decomposition of the one-loop four-gluon amplitude is 



E N c Tr (TMU T a <» T a °^ T a ° (4) )^JJi (<r (1) , <r (2) , a{3) , <r (4) ) 

+ E Tr(T a -( 1 )r a -( 2 ))Tr(r a -( 3 )r a -( 4 ))A4 ;3 (o-(l),(T(2);o-(3),(T(4)) . 
<xeS 4 /zf 

+ n } Yl Tr(r a -(i)T a -( 2 )T a -(3)r aCT ( 4 ))^{ 2l (cj(l),cj(2),cj(3),cj(4)) 
0-GS4/Z4 

+ n s E Tr(T a -( 1 )TM2)T a -(3)T a -( 4 ))A! t ° 1 1 (cj(l),(T(2),cj(3),a(4)) 

(2.5) 

In the first term, the permutation a lies in the set of all permutations (£4) of four objects with 
purely cyclic ones (Z4) removed. In the second term, a is again in the set of permutations of four 
objects but with two factors of Z2 removed, corresponding to exchanging the indices within each 
trace, as well as another Z2 removed corresponding to interchanging the two traces. The leading 
partial amplitude is indicated by the subscript '4;1', and the sub- leading partial amplitude 
is indicated by the subscript £ 4;3'. In this process, only a single sub-leading color structure 
appears. Processes with more external legs can have more sub-leading color structures. Note 
that in eq. (^j|), we have abbreviated the dependence of the A n -j on momentum k\ and helicity 
A; by writing the label I alone. 

The superscripts [1], [1/2] and [0] label the spin of the particles circulating in the loops. 
The third and fourth terms are part of the leading partial amplitude and correspond to the 
contributions of nj flavors of fundamental representation quarks and n s flavors of complex 
scalars in the (iV c + iV c ) representation. In QCD, below the top threshold rif = 5 and n s = 0. 
Although fundamental scalars do not appear in QCD, it is convenient to keep explicit dependence 
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on the number of scalars; when the number of bosonic states are equal to the number of fermion 



states certain supersymmetric Ward identities [39] must be respected which can be used as a 
check on results. 

In our convention, each matter representation has 4N C color and helicity degrees of freedom. 
This does not correspond to the standard assignment of 2N C degrees of freedom for scalars, but 
has the advantage of making supersymmetry identities more apparent since the number of states 
matches that of a Dirac fermion. 

Similarly the five-point amplitude is, 



^l-loop 5 



£ N c Tr(a(l)...a(5))4l(a(l),...,a(5)) 

+ Yl ^Kl)a(2))Tr( f 7(3) ( r(4) f 7(5))A 5; 3(a(l), <r(2); ^(3), <r(4), a(5)) 
°es 5 /(z 2 xz 3 ) 

+ «f E TV(a(l)...a(5))4^ 21 Kl),...,a(5)) 



(2.6) 



+ n s Yl Tr(a(l)...a(5))Afj l (a(l),...,a(5)) 

<r&S 5 /Z 5 

where S5/Z5 is the set of non-cyclic permutations of five objects and S^jiZ^ X Z3) is the set 
of permutations of five objects that do not leave the product of two and three element traces 
invariant. 

In the pure external gluon case the leading color amplitudes (broken down by the spin of 
the loop particle) play the role of primitive amplitudes in the sense that all sub-leading partial 
amplitudes can be obtained by the appropriate permutation sum given by 



A 



n;c>l 



(1,2, 



l:c, c + 1, 



-1 



,c-l 



(2.7) 



a£COP{a}{/3} 



where ctj G {a} = {c — 1, c — 2, . . . , 2, 1}, $ € {/?} = {c, c+ 1, . . . , n — 1, n}, and COP{a}{f3} is 
the set of all permutations of {1, 2, . . . , n} with n held fixed that preserve the cyclic ordering of 
the oti within {a} and of the (3i within {/?}, while allowing for all possible relative orderings of 
the a.i with respect to the Pi. For example if {a} = {2, 1} and = {3, 4, 5}, then COP{a}{(3} 
contains the twelve elements 



(2,1,3,4,5), (2,3,1,4,5), (2,3,4,1,5), (3,2,1,4,5), (3,2,4,1,5), (3,4,2,1,5), 
(1,2,3,4,5), (1,3,2,4,5), (1,3,4,2,5), (3,1,2,4,5), (3,1,4,2,5), (3,4,1,2,5). 



(2.8) 



The color decomposition of the qqgg one- loop amplitude is [24, 40 1 



N c ]T (TM3)rM4)) i2 H A 4; i(l ? -,2 ? ; ( r(3),a(4)) 

<tGS 2 



+ Tr (T a3 T a4 ) 5 J 1 A 4 ; 3 (l $ ,2 9 ;3,4) 



(2.9) 
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Similarly, the color decomposition of the qqggg amplitude is [19, 24] 



A(l 5 ,2 g ,3,4,5) = g 5 N c £ (T<M3) nw^K^ 1 i 5; i(lg, 2 g ; <t(3), a(A), a(S)) 

+ J] Tr(T<M3)T<M4)) (T^W) * A 5 . 3 (l ?! 2 g ; a(3), <x(4); a(5)) (2 1Q) 

+ £ Tr(TM3) T a 4rM 5 )) ^ ^. 4 (i fj 2g ; a(3), 4, <r(5)) 

tTG5 2 

In the partial amplitude an additional semicolon separates the gluon sandwiched between 
the quark indices (the last gluon in A^) from the other two gluons. 

For the case of external fermions the leading partial amplitudes cannot play the role of 
primitive amplitudes. The primitive amplitudes with nearest neighboring quarks are given by 
subdividing A n -\ into smaller pieces pf], 



An;l(lq, 2q] 3, . . . , n) — A n ' ^ ^ (lq, 2q, 3, . . . , Tl) -^n ' ^ ^ 2gj • • • j 



n 



, H£ A L,ll/2](,_ n o 
N 



(2.11) 



J(l ? ,2„3,...,n) 



All color factors have been extracted into the coefficients of the primitive amplitudes. The 
primitive amplitudes themselves are independent of the numbers of colors. The labels L and 
R on the primitive amplitudes refer to whether the fermion line which enters a diagram turns 
either 'left' or 'right'. For the purposes in this paper we may take eq. (2.11) to be the defining 
equation for the two-quark n — 2 gluon amplitudes, when all gluons are between quark and the 
anti-quark in the cyclic ordering. Further details, and the definition of the primitive amplitudes 
when gluons are also between the anti-quark and quark, may be found in ref. |f24f| . The L and 
R primitive amplitudes are related by an inversion of the ordering of legs, 

A^ J \l q ,2,..., k -l,k q ,k + l,...,n) = {-l) n A^ J \l q ,n,...,k + l,k q ,k-l,... ,2) . (2.12) 

For the L, [1] primitive amplitudes all gluon legs between the anti-quark and the quark in the 
cyclic ordering are attached to the fermionic part of the loop, while all gluon legs between the 
quark and the anti-quark are attached to the bosonic part of the loop. (See ref. [^4|] for further 
details.) 

The sub-leading partial amplitudes of two quark processes are given by permutation sums 
similar to eq. fl2.7|) for all-gluon processes, 



Ai;oi(lq, 2 q \ 3, . . . , c+1; c + 2, . . . , n) = 

(-ir 1 E 

aeCOP{a}{/3} 



(2.13) 



where a; G {a} = {c + 1, c, . . . , 4, 3}, G {/?} = {l q , 2 q , c + 2, c + 3, . . . , n — 1, n}, 
For example, for the four-point qqgg amplitude, 

A 4 . 3 (l q , 2 q ; 3, 4) = Aj' [1] (l q , 2 q , 3, 4) + a£" [1] (l q , 2„ 4, 3) + A±' [1] 3, 2 q , 4) 
+4' [1] (1 5 , 3, 4, 2,) + ^' [1] (l q , 4, 3, 2 9 ) + 4' [1] 4, 2 q , 3) , 



(2.14) 
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Symmetry relations among nj and n s terms cause them to cancel out. 
The five-point relations are, of course, a bit more complicated, 

A 5;3 (l g -,2 g ;4,5;3) = £ Aj" [1] 2„ a(3), <r(4), <r(5)) + 4' W 0-q, 4, 2„ 5, 3) 

+4' [1] 4, 2 q , 3, 5) + 4' [1] (lg, 5, 2„ 3, 4) + A%> [1] 5, 2 q , 4, 3) ( 2 ' 15 ) 
+A L b [1] 4, 5, 2„ 3) + 4' [1] 5, 4, 2„ 3) , 



45. 4 (1 ? ,2,;3,4,5) = E 



-A*" [1] 2, , a(5), (7(4), a(3)) - A^' [1] a(5), 2„ a(4), <r(3)) 



- Af (l f , a(5), a(4), 2, , a(3)) - Af (1 ?J a(5), (7(4), a(3), 2,) 

-S- (4' 11/21 (Iff, 2 ?J (7(3), (7(4), (7(5)) + [1/2] (1 5 , (7(3), 2 ? , (7(4), (7 



(Aj- M (l ff , 2„ a(3), a(4), a(5)) + Ag' M (l ff> a(3), 2„ a(4), er(5))] 

(2.16) 

The primitive decomposition of four quark amplitudes have not been explicitly given in 
the literature. Nevertheless, such a decomposition can be performed by following the ideas in 



ref. [24|. The first step in performing a decomposition into primitive amplitudes is to perform 
the usual decomposition into coefficients of independent color factors. Then as in the previous 
cases one defines a set of primitive amplitudes in terms of a set of 'parent' diagrams. As for the 
other cases, the essential distinguishing feature of the primitive amplitudes is that the colored 
external legs of all contributing diagrams have a fixed cyclic ordering. It is this general feature 
of primitive amplitudes which we use in subsequent sections to present a systematic description 
of the soft properties of the one-loop amplitudes. 



3 Review of collinear and soft properties of QCD amplitudes 

The properties of QCD tree amplitudes in the limits where momenta become collinear or soft has 
been extensively discussed in the literature. (See for example ref. |0].) Similarly, the properties of 
one-loop amplitudes in collinear limits have also been extensively discussed |2|, |||, |26|, through 
O(e ). Although, the properties of one- loop QCD amplitudes as external momenta become soft 
has been less extensively discussed, it is straightforward to extract explicit expressions for the 



behavior using the known four-point 41] and five-point amplitudes [22, 23, ^J]. In this paper 
we wish to extend these results to higher orders in the dimensional regularization parameter, so 
that they can be inserted into NNLO computations. 

In addition to the application of using the soft and collinear splitting amplitudes to evaluate 
the infrared divergent regions of phase space integrals, the understanding of the collinear and 
soft properties of one-loop amplitudes has led to a method for constructing amplitudes from 
their analytic properties. This has been used to simplify the computation of the Z — > 4 parton 
helicity amplitudes and to obtain infinite sequences of one-loop maximally helicity violating 
amplitudes |, ||, ||. 
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We now review the known soft and collinear properties of the one-loop amplitudes in prepa- 
ration for our extension of the results to all orders of the dimensional regularization parameter. 



3.1 The collinear behavior 

As the momenta of two massless particles become collinear, one-loop amplitudes factorize in 
the way depicted in fig. HI. As discussed in ref. [33], for infrared divergent massless theories the 



situation is more subtle than in infrared finite cases since the various contributions to individual 
integral functions diagrams may not have a smooth behavior as the intermediate momentum 
becomes massless. In particular, in terms of the diagrams of covariant gauges, the picture of 
fig. |l] is slightly misleading since there can be contributions from diagrams with no explicit 
propagator on which to factorize the amplitude; the kinematic poles arise instead from the loop 
integral. (This subtlety is related to the interchange of the limit e — > with that of a kinematic 
variable vanishing s a b — > 0.) Nevertheless, the complete amplitudes do obey the factorization 
described by fig. [|. 

Figure 1: The generic behavior of one- loop amplitudes in the limit as two external momenta become 
collinear. The shaded disc represents the sum over tree diagrams and the annulus the sum over one-loop 
diagrams. 

More explicitly, as two adjacent momenta k\ and k 2 become collinear, we may factorize an 
n-point tree amplitude on the s%2 kinematic pole in terms of a three vertex and an (n — l)-point 
amplitude, 

AT (T2, . . (^E ^^f^ ) di(K) AT -^ KX ^ ■ ■ ' (3 ' 1} 

where A specifies the polarization and k\ = zK and k 2 = (1 — z)K with K = k\ + k 2 . This basic 
structure is independent of the particle type, although we have written eq. ([O]) for the case 
of an intermediate gluon. The values of the p tree functions do however depend on the particle 
types and are simply given by the three-vertices 

^4T 92 = »V2(fcfr • e 2 + k 2 ■ e x e» - h ■ £2 e£) , 

^-.u, tree * _ a 



-TV/', tree 



In the first two cases \x represents the Lorentz index of the intermediate gluon which in the last 
case is replaced by a spinor index. Although it makes no difference in the results, we use the 



non- linear Gervais-Neveu [43] gauge since the vertices are particularly simple. 

After inserting an explicit representation of the helicity states [pj], the collinear limits of 
tree amplitudes may be re-expressed as (see e.g. ref. 

< rc °(l Al ,2 x \...,n)^ ]T Split toc A c (l Al , 2 A2 ) (K\ 3, . . . , n) . (3.3) 

A=± 
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Figure 2: The three vertices describing (a) £>£l^f S2 , (b) I>£i^° a and (c) 2>£* r £ f 



The splitting amplitudes in eq. ( p.3| ) have square-root singularities in the collinear limit. For 
convenience we have collected the tree-level helicity splitting amplitudes in appendix |A]. For most 
computational purposes, the helicity form in eq. (3^3) is the more convenient one, but if one is 
working entirely in the context of conventional dimensional regularization where the polarization 
vectors become (4 — 2e)-dimensional, the representation in terms of formal polarization vectors 
and spinors is also useful. 

The behavior of one-loop primitive amplitudes as the momenta of two adjacent legs becomes 
collinear is similar and is given by 



4r ioop (i, 2, . . . , n) ±H 2>&~ U £ o^r/n^ r p (K\ 3, . . 



A 

-A/ T^\^\l 



^fi, 1-loop / • \ - £ M ( K ) £ u( K ) \ 9 trcc A 

In the spinor helicity language, this becomes [26. 



(3.4) 



^i-ioop^Aj ^ X \..., n )^ {split^(l Al , 2 Aa ) A£?F(K\ 3, . . .) 

+ Split 1 ;} 0015 (l Al , 2 Aa ) A^\(K x , 3, . . .) 



A=± 



(3.5) 



where A represents the helicity, A^ n ^ v and A^ ee are one-loop and tree sub-amplitudes with a 
fixed ordering of legs with legs 1 and 2 consecutive in the ordering. 

In section ^| we give expressions for the one-loop splitting functions P^' 1-1000 and Split 1 " 1000 
to all orders in the dimensional regularization parameter. 



3.2 Soft behavior 

In the limit that a gluon momentum becomes soft we may again factorize tree amplitudes on 
kinematic poles (see e.g. ref. In taking k± — ► the two kinematic poles in the color ordered 
tree amplitude which may diverge are l/s±2 and l/s n i since n, 1 and 2 are consecutive legs in 
the cyclic color ordering. This yields the factorization relation 

< ee (l ± , 2, . . . , n) ^V ree (n, 1* 2) A^(2, 3, . . . , n) . (3.6) 

In terms of polarization vectors, the tree-level soft function, 5 tree , is the familiar eikonal factor, 



2ef ■ k n 2et ■ k 2 



Snl 



S12 



(3.7) 
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In spinor helicity notation, the soft amplitudes are, 



-.tree 



(n,l+ 2) 



<"2) 
(nl)(12) 



S tree (n,l~,2) 



[nl] [12] 



(3. 



The soft amplitudes are independent of the helicities and particle types of the neighboring legs 
n and 2. 

The behavior of one-loop gluon amplitudes in the soft limit is similar to the collinear behavior. 
As the momentum k\ becomes soft the primitive one-loop amplitudes become 

^- lo °P(l, 2, . . . , n ) k ^°S tTee (n, 1+2) A^°° p (2, 3, . . . , n) + S lAoop {n, 1 ± , 2) A^_\(2, 3, . . . , n) . 

(3.9) 

The onedoop soft function may be extracted from four- 4C, 41, and five-parton |^, [23|, 
one-loop amplitudes that are known through O(e ), with the result, 



S 



l-loop 



71,1^2) 



r 1 



H 2 {-Sn2) 

(s n l)(-si2) 



+ 



7T- 

6~ 



(3.10) 



As for the collinear case, in section we will generalize this to all orders in the dimensional 
regularization parameter. We will also generalize the discussion to include the case of one-loop 
primitive amplitudes with external fermions and present a representation in terms of formal 
polarization vectors. As we shall see, although the soft function does not depend on the neigh- 
boring external particles in the primitive ordering of legs, it does depend on whether the soft 
gluon is attached to a gluonic or fermionic part of the loop. 



4 One-loop soft and collinear splitting amplitudes 

In this section, we will apply the analysis of collinear limits and their relationship to infrared 



divergences of ref. [33 1 to obtain explicit expressions for the splitting amplitudes to all orders in 
the dimensional regulator, e. We shall also present explicit expressions for the soft amplitudes 
to all orders in e. 

In our organization of the results, there are two contributions to the one-loop splitting 
amplitudes; the 'factorizing' and the 'non-factorizing' pieces, Split 1 " 10013 = Split fact + Split non " fact . 
The factorizing contributions are those which one naively expects for one-loop splitting. That 
is, they are derived from those diagrams, shown in fig. ||, in which one particle from a tree-level 
amplitude splits into two via a loop. They are called factorizing contributions because they 
factorize on a single particle propagator; one can amputate the splitting term from the diagram 
by cutting a single line. 

The non-factorizing contributions are quite different in nature. These arise from non-smooth 
behavior of infrared divergent integrals as a kinematic invariant vanishes. This non-smooth 
behavior arises because we take e — > before we take the kinematic invariant to vanish. These 
non-factorizing contributions are proportional to the tree-level ones since the infrared singular 
parts of one-loop amplitudes are proportional to the tree-level amplitudes. 

Although the division of the splitting amplitudes into factorizing and non-factorizing pieces 
contains some arbitrariness and is gauge dependent, the sum of the contributions is gauge 



10 



+ 1 





+ x> + 




Figure 3: The diagrams in a massless theory (ignoring tadpoles) that need to be calculated to obtain 
the factorizing contribution to the loop splitting function. The dotted line represents the off-shell leg on 
which the collinear factorization is performed. 

independent. The calculation of the diagrams in fig. |3] always gives the complete non-singular 
contribution to the one-loop splitting amplitudes, but, depending on the gauge choice, also gives 
some portion of the infrared singular contribution. However, since the divergence structure of 
one- loop splitting amplitudes is known (see appendix [B| for their values through C(e )) and since 
each 1/e pole is uniquely identified with a known 'discontinuity function' (see reference |33| and 
section below) we can push all singularities into the non-factorizing piece, and define the 
factorizing piece to be the infrared finite part of the diagrams in fig. ||, not involving any of the 
discontinuity functions. The non-factorizing contribution may then be determined by finding 
the linear combination of discontinuity functions that matches the known infrared divergence 
structure of the amplitudes. 

We will first compute the diagrams that yield the factorizing parts of the one-loop splitting 
amplitudes. We then use the known structure of the infrared singularities of the leading partial 
amplitudes to derive the non-factorizing contributions to one-loop splitting. By combining 
the factorizing and non-factorizing parts we obtain the full one-loop splitting amplitudes for 
the leading partial amplitudes. Finally, we will derive the soft factorization properties of one- 



loop amplitudes. Our results will be in terms of the bare splitting amplitudes; in section 
we give the appropriate ultraviolet subtraction for obtaining the splitting amplitudes for the 
renormalized amplitudes. 



4.1 Factorizing contributions to the one- loop splitting amplitudes 

The factorizing contributions to one-loop splitting amplitudes are determined by computing 
triangle and bubble diagrams like those shown in fig. 0. Any poles in e that appear must come 



from one of two discontinuity functions [33] 

/i(si2,e) 
hi *.>■<) 



H 2 



e(l-2e) \-s 12 

i ( » 2 V 



(4.1) 



-S12 



After removing the poles by subtracting off terms proportional to these discontinuity functions 
the factorizing contribution to the splitting function is determined by inserting a complete set 
of states on the (off-shell) fused leg, imposing the collinear limit and then taking the fused leg 
on-shell. 
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The results below are computed in the Feynman background field gauge and apply to an 
SU(N C ) gauge theory with matter content of rif flavors of massless Dirac fermions in the funda- 
mental representation and n s flavors of massless complex scalars in the (N c + N c ) representation. 
In each case, we describe a parton splitting into two final state partons labeled 1 and 2 which 
carry momentum fractions z and 1 — z respectively. 

The difference between dimensional regularization schemes is parameterized by the quantity 
Sr. In the dimensional reduction [46| or four-dimensional helicity [^TJ schemes used to compute 
one-loop helicity amplitudes, both external gluons and those inside loops, have two polarization 
states. In these schemes, 5r = 0. In the 't Hooft-Veltman scheme, external gluons have 
two polarization states, but those inside loops have [e]-helicities, i.e. polarization states that 
point into the extra — 2e dimensions. For this scheme, 5r = 1. In conventional dimensional 
regularization all gluons, external and internal, have [e]-helicities. The [e]-helicities of the internal 
gluons are accounted for by setting 8r = 1; those of the external gluons must be accounted for 
by spin sums over .D-dimensional polarization vectors. 



4.1.1 Factorizing contributions to one- loop g — > gg splitting amplitudes 

The result of the triangle and bubble graphs is (after stripping the coupling and color factors) fl33| 



f)[t, fact 
9^9192 



Tit I'-^-i + iJ^-^ 



(ei-£2) 



(ki-k 2 



where 



and 



6 



-S12 



ev- 



er 



(l-2e)(2-2e)(3-2e) ' 

1 r(i + e )r 2 (i-e) 



(47T) 



2-e 



rfl-2e) 



(4.2) 



(4.3) 



(4.4) 



In background field Feynman gauge, which is a convenient gauge for evaluating the diagrams, 
^9^9192 ^ s fi n it e due to a Ward identity, so no discontinuity functions need to be subtracted off. 
The factorizing part of the one loop splitting function, Split fa ^ t (l Al , 2 A2 ) (where the A's label 
helicities and are defined as if all particles exit the diagram), is derived from D by 



i fact 
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•r-\^, fact 
9^9192 



(Ai,A 2 



ejf{K)el{K) 



2ki ■ k 2 

A 



^tree / \ 

dei(K) An - l[ --- K 



(4.5) 



Split fa ^ t (l Al , 2 Aa X r f c 1 (. . . K 



This yields 



TV 



n 



Splie?(l^,2^) = -^[i-e5 R -^ + jf- 



2V2(k r k 2 ) 



(ei-^2) 



(ki-e 2 )(k 2 -ei) 



{ki-k 2 



(4.6) 
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where K = k\ + hi- Expressed in terms of spinor inner products 44 , |]] the result is that 



Split 



Split 
Split!?** 



6^/1-0-) 



3 



iV c A c 



z(l - z) 



3 V" "~" iV c N c 
,2") = Split^ ct (l-,2+) = 0, 



(12) 



(4.7) 



2] 



with the remaining terms, Split_(l Al , 2 A2 ), given by parity inversion. Except for Split + (1 + 
and Split_(l~, 2 _ ) which vanish at tree level, the factorizing contributions to the one 
splitting amplitudes are proportional to the tree-level splitting amplitudes, 



Split f ^(l Al ,2 A2 ; 



crSplit^ 6 (l Al , 2 Aa )r|"" 9192 ' fact (-A, l Al , 2 



X) 

loop 
(4.8) 



where 



r ™?2,fact (A)1±)2T)=()) 



„9^3i92,fact 
S 1 



±,1 T ,2^) 



Sl2 



l-eS R 



nf n s \ 2z(l — z) 

N~ C + N~J (l-2e)(2-2e)(3-2e) 



(4- 



4.1.2 Factorizing contributions to one-loop g — > qq splitting amplitudes 

For g -^qq splitting, the sum of the color and coupling stripped triangle and bubble graphs is 

' .2 \ e 



fi, fact 



13 



Sl2 



+ 



1 



+ 



(1 - Sr) 



6e(l-2e) 3(1 - 2e)(3 - 2e) (1 - 2e)(2 - 2e)(3 - 2e) 



1 

i 

~N~ n 



1 



e 2 2e(l 

2n/ + n s 



+ 



+ 



2e) 1 - 2e (2 - 2e) 
2(n/ - n s ) 



(4.10) 



3e(l-2e) 3(1 - 2e)(3 - 2e)_ 

Here we see a number of terms that are proportional to the discontinuity functions f\ and f%. 
After removing these terms, we have, 

■ 2 V r 5-2e-3fo 



2> 



fact 



'V2 



-S12 



1 l + <5 fl (l-2e) 



3(1 - 2e)(2 - 2e)(3 - 2e) 
rif — n s 2 



A c 2 (l-2e)(2-2e) A c 3(1 - 2e)(3 - 2e) 



(4.11) 



After inserting a complete set of helicity states in an expression analogous to eq. (| 
find that the one loop splitting function is proportional to the tree-level result for each 

as in eq 

,2 



WC 

allowed 



helicity configuration. With r ^ <?1<?2 ' fact defined 



n 9^gil2, fact 



(A, 1 , 2 T ) 



Sl2 



5 - 2e - 35 R 



, we have, 

1 l + 5 R (l-2e) 



+ 



3(1 - 2e)(2 - 2e)(3 - 2e) A" 2 (1 - 2e)(2 - 2e) 
nf — n s 2 



+ 



3(1 -2e)(3-2e) 



(4.12) 
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4.1.3 Factorizing contributions to one-loop q —* qg splitting amplitudes 



For g — > qg splitting, the sum of the coupling and color stripped triangle and bubble graphs is 



V 



, fact 



] V2 



Sl2 



unh(ki-e 2 ) 
(fci-fca) 

iV 2 + 1 1 - 



iV 2 e(l-2e) iV 2 (l-2e)(2-2e) 



After subtracting out terms proportional to the discontinuity functions, we get 



2? 



i, fact 

q^qigi 



V2 



"Sl2 



^lfe(fcr£2) 



N? + l 



1-6(5 



iV 2 (l-2e)(2-2e) 



Following the same procedure as above, we find that 

r ™, 2 ,fact (±)lT;2±) =Qj 



-312 



1 ~t~ "T777 



1 \ (l-z)(l-e6 R ) 



N*J (l-2e)(2-2e) ' 



(4.13) 



(4.14) 



(4.15) 



For g — ► g\q 2 splitting, we need simply interchange z and 1 — z. 



4.2 Non-factorizing contributions 

In ref. |33| it was shown that all non-factorizing contributions may be linked to infrared diver- 
gences which have a universal structure for an arbitrary number of external legs. The coefficients 
of the infrared divergences may be used to fix the coefficients of all integral functions associated 
with non-factorizing contributions to the splitting amplitudes. 

The integral functions to be used are selected by demanding that they contain no infrared 
divergences other than ones which may appear in the splitting amplitudes. The infrared di- 
vergences of the integral functions may be obtained from the explicit forms of the integrals 
contained in, for example, refs. [^, |4^| . By systematically stepping through the list of all inte- 
grals and their discontinuities one can construct a list of infrared divergent functions (and their 
associated higher order in e parts) that may appear in the soft or collinear splitting amplitudes. 

More explicitly, the complete one-loop splitting amplitudes are given by 

Split 1 - l0 °P(l, 2) = Spht fact (l, 2) + ci/i(*i 2] e) + c 2 f 2 (s 12 , e) + c 3 / 3 (s 12 , e, z) + c 4 / 4 (si2, e, z) , 

(4.16) 

where the Cj are coefficients to be fixed using known infrared divergences and the discontinuity 
functions are (see appendix |C]) 

h(s 12 ,e) = -i^Msn) = TT^i—T 

e(l - 2e) v -si2 / 

/2(si2,e) = -i/U 2e J 3 (si2) = -%{-^— ) , 
h(si2,e,z) = -/20l2,e) + lim ^ 2e s n is 12 ll^ 

fel ||«2 

fA(si2,e,z) = -f 2 (s l2 ,e) + lim ^/i 2e s 12 s 23 ll™ 

fel ||«2 
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m=l 
oo 

E e m Li„ 



m=l 



1 - Z 

—z 

—z 
l-z 



(4.17) 
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where Li m (x) is the mth polylogarithm 



Lii(x) 
Li m (x) 



ln(l - x) 

x dt 

— >->l m -i(t) 

o t 



(m = 2,3,...) 



£ 

71=1 



(4.18) 



The scalar integrals appearing in the discontinuity functions are 



(a) J 2 (si 2 ) = 

(b) X 3 (si2) = 

(c) T\%{s nl ,s 12 ) 

(d) ll%(s 12 ,s 23 ) 



d 



4-2e„ 



1 



( 27r )4-2, p 2( p _ fcl )2( p _ fcl _ fe2 )2 ' 
d 4 " 2e p 1 



(4.19) 



( 27r )4-2e p 2 (;p _ kl) 2( p _ h _ fayfr + K y 



d 4 ~ 2t p 



1 



(2vr) 4 - 2e p 2 (p - ki) 2 (p — k\ — k 2 ) 2 {p -k 1 -k 2 - k 3 ) 2 ' 

These integrals are depicted in fig. ||. No other integrals appear because they are either 
smooth functions in the collinear limits, or they do not have the correct infrared divergences to 
match those appearing in the splitting amplitudes |33[]. It is easy to see from eq. ( [1.17 ) that 
/4(si2> £) z) = f 3 {s\ 2 , e, 1 — z). This relation is apparent even in the integral function by shifting 
the loop momentum of T\^{s\ 2 , s 23 ) such that p — ► — p + k\+ k 2 . 

Xzx 2 

(a) 1 




n- 1 



(c) 



77 



(d) 



Figure 4: The scalar integrals that can appear in the collinear splitting amplitudes. Legs 1 and 2 are 
collinear. 



To fix the coefficients c% we proceed as follows. An n-point primitive amplitude will have a 
known singular structure of the form [18, ||, |i~9| j, 



A 1-loop 



singular 



tree 



+ 



(4.20) 



where the coefficients Sj and depend on the particular n-point amplitude under consider- 
ation. Consider now the collinear or soft limits given in eqs. ( [3,5| ) and (p.Sj). The n — 1 point 
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amplitudes appearing on the right-hand-side of these equations can have a different set of singu- 
larities. The difference between these two sets of known singularities must then be absorbed into 
the splitting or soft amplitudes Split 1_loop or s 1 ' 10 ^. More directly, we can use the known results 
for the splitting amplitudes through O(e ), which we have collected in appendix [B], to determine 
the singularities. Using table p] we can then uniquely determine all Cj coefficients appearing in 
eq. ( [4.16 ). The unique 1/e singularities of the functions are listed in the first column of the 



table. The second column lists the discontinuity function associated with each singularity. One 
adjusts the q coefficients so that the divergences in the splitting amplitudes are correct. This 
uniquely fixes the splitting amplitudes to all orders in the dimensional regularization parameter 

6. 



Singularity 


Non-Factorizing Contribution 


cr 

e 


/i(si2,e) 


crln(-s it i +1 ) 
e 


/2(si2,e) 


cr ln(z) 
e 


/30l2,e,2) 


cr ln(l — z) 
e 


/ 3 (s 12 ,e, 1 - z) 



Table 1: The potential 'non-factorizing' contributions to the splitting amplitudes in the k\ \\ ki channel. 



Since the divergences of the primitive amplitudes are always proportional to tree amplitudes, 
the non-factorizing contributions to one-loop splitting amplitudes will all be proportional to the 



tree-level splitting amplitudes. In analogy to eq. (4.8), we write 



Splie° A n - fact (l Al , 2 Aa ) = c r Split^ (l Al , 2 A2 ) r g°n-fact 



(4.21) 



for each type of splitting function. There is no need to specify the helicity structure of r<j on ~ fact 
because the factor is the same for all allowed helicity configurations. For the T> functions we 
can write 



V 



a. non-fact .tn/i. mv w< 

£-1,2 = c ? V K->l,2 r S 



fi, tree non-fact 



(4.22) 



4.2.1 Non-factorizing contributions to one loop g — gg splitting amplitudes 

In order to demonstrate the procedure more explicitly, we extract the non-factorizing contri- 
butions to one-loop g — gg splitting using pure gluon amplitudes. The unrenormalized pole 
structure of the leading n gluon partial amplitude A n: i(l, 2, . . . , n) is 



singular 



-cr^ 9 ' tree 



1 n 
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3iV c 



1 n s 
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(4.23) 
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where we identify the label n + 1 with 1. In the limit that gluons 1 and 2 become collinear, we 
verify that for the helicity configurations forbidden at tree-level, 



Splits-^- non - fact (±, 1± 2 ± ) = , 
while for all configurations allowed at tree-level 



„9->9l92, non-fact 
r S 



1_ ln(-g 12 ) ln(£) ln(l - 
e 2 e e e 



+ O(e ) 



(4.24) 



(4.25) 



This may also be read off from the known results [26, 24 1 for the splitting amplitudes through 
O(e ) which we have collected in appendix [B|. Focusing on the 1/e poles we use table [l] to deduce 
that 

r ^ 5l52 , non-fact = f ^ g) + f ^ ^ z) + ^ 1 _ ^ ^ ^ 

where the discontinuity functions /j are defined in eq. ( [4.17 ). The combination of functions 
appearing here may be further simplified using eq. ( |C.10|) . 



4.2.2 Non-factorizing contributions to one-loop g — > qq splitting amplitudes 



To extract the non-factorizing contributions to one-loop g — > qq splitting we need the unrenor- 
malized pole structure of the leading all-gluon partial amplitudes, eq. ( 4.23 ), as well as that of 
the leading two quark partial amplitudes A n: i(lq, 2 q , 3, . . . , n), 



A 2q 



singular 



-c r 
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iVf e 2 



-«12 
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(4.27) 



From this it follows that 

ff — *9l92; non-fact 

11 1 A, v , v ln(-si 2 ) 13 3 



2n/ 
3 



+ O(e ) . 



(4.28) 



As before, this may also be obtained from the known results for the splitting amplitudes through 
0(e°) which we have collected in appendix [B[ Again, focusing on the 1/e poles, we read off 



r.g^ll 92, non- fact 



/3(«i2, e, z) + h{si2, e, 1 - z) - -tto h{svz, e 
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iV 2 ' 

c 

13 3 2n/ n s 
¥ + 2ivf ~ 3iVc ~ 3iVc 



(4.29) 



/i(si2,e) . 



4.2.3 Non-factorizing contributions to one loop q^> qg splitting amplitudes 



We need only the pole structure of the leading two quark partial amplitudes, eq. ( |4.27| ), to derive 
the non-factorizing contributions to one-loop q — > qg splitting. From 

„™ 2 ,non-fact = _ JL + I (hi(- Sl2 ) + ln(l - z) - ^) + 0(e°) , (4.30) 



e 2 e 



we obtain 



„™ 2l non-fact = g) + ^ j _ ^ _ _L/ g ( ai2j e , z ) 



(4.31) 



Again, for g — * g\qi splitting we simply interchange z and 1 — z. 
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4.3 Full one-loop splitting amplitudes for leading color partial amplitudes to 
all orders in e. 



We now combine the factorizing and non-factorizing terms to give the full one loop splitting 
amplitudes for leading partial amplitudes to all orders in e. 

For g — ► g\g 2 splitting, there is one independent helicity configuration that is non-zero at 
one loop but vanishes at tree-level, 

Split™' 1 - loop (l+,2+) = 



V V ' (12? 1 



2c r 



(12) 2 (l-2e)(2-2e)(3-2e) \-s 12 ) 
We present the rest of the splitting amplitudes in terms of rs- 
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In terms of formal polarization vectors, 
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(4.34) 



For g ^> qiq 2 splitting, we obtain the same result for each allowed helicity configuration 
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(4.35) 

For g — > ^1^2 splitting we simply interchange z and 1 — z. In fact, the 5 — > and g ^ gg 
splitting functions are symmetric in z — > 1 — z, as can be seen from the results of appendix |C|. 
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In terms of formal polarization vectors the T> function is 



2#5& P = ^u2^v 1 c r rr' iq2 ' 1 - loop (x,i ± ^: 



(4.36) 



Finally, for q — ► qig 2 and q — > g±q 2 splitting, 
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(4.37) 

For q — * g\q 2 and q — * q\g 2 splitting we simply interchange z and 1 — z. In terms of formal 
polarization vectors, these are written as 



1-loop _ ' 
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(4.38) 



4.4 Renormalizing the one-loop splitting amplitudes 

We can easily convert to the renormalized splitting amplitudes by performing an MS ultraviolet 
subtraction^] using the fact that the splitting amplitudes implicitly contain one power of the 
couplin J] and that the ultraviolet divergence of an n-point one-loop amplitude is given by 



A 



l-loop 



,tree™- 2 f 11 1 n f 1 n s 



= Cr ^ ree 7-;?-5F ■ ( 4 - 39 ) 

UV divergence e \ O 3 JS c D iV c 



All other poles in e in A\ loop are infrared singularities. Since all of the singular terms have been 
pushed into the non-factorizing term, renormalization can be accomplished by altering just that 

^The MS subtraction jl9j is precisely defined only to O(e ). Beyond O(e ) there are several different definitions 
pl| , pK| pl| and [^(J (which contains two different definitions in sect. 4 and 5). We choose the second definition of 

^The MS running coupling depends on the regularization scheme in which the one- loop amplitude is renormal- 



ized; e.g. in the 't Hooft-Veltman scheme, renormalizing the amplitude through eq. (1.40) amounts to replacing 
the bare coupling a s with the running coupling af v . The difference between running couplings in the different 
schemes is of 0(e°) §(], |§. 



19 



piece: 

a— >6iC2, non-fact a— >feiC2, non-fact 1 / 1 ' 1 ^/ 1 \ i . 

Ts >™ ~ rs ~ ~e [~6 " 3N~c " 6N~J ■ (4 ' 4Uj 

Following this procedure we can rewrite any of the splitting or soft amplitudes presented in this 
paper so that they hold for renormalized instead of bare amplitudes. 

4.5 Primitive decomposition of one-loop splitting amplitudes 

As discussed in section 0, one-loop amplitudes are naturally broken down into color ordered 
partial amplitudes, which can themselves be broken down into gauge invariant primitive com- 
ponents. Once we determine the splitting rules for the primitive amplitudes, then formulae like 
eq. ( p.7D can be used to generate the splitting rules for any partial amplitude. 

In fact, collinear splitting is a property of full partial amplitudes, not just of the primitive 
ones. Furthermore, it is only the leading partial amplitude, A^-i which has the same color struc- 
ture as the tree-level amplitude, that receives splitting contributions from both Split tree and 
Split 1 " 10013 . The sub-leading partial amplitudes, which have no tree-level counterparts, require 
only Split tree . The primitive components of the sub- leading partial amplitudes do receive split- 
ting contributions from Split 1 " 10015 , but these contributions cancel in the permutation sums that 
build the sub-leading partial amplitudes. This can easily be seen for the all-gluon case using 
eq. @, 

A n .j(l, 2, . . . , j - 1; j, . . . , n) ±E (-1)^ 1 Split trcc (l, 2) £ A% ] _ 

*ecop{a}{p} ^ 41 , 

+(-ir i spht i - io °p(i,2) y, A n-ii°), 

aeCOP{a}{!3} 

where oti 6 {a} = {j— 1, j— 2, . . . , 3, c}, £ {f3} = {j, . . . , n— 1, n}, and c is the coalescence 
of gluons 1 and 2. The sum of tree-level amplitudes over COP permutations vanishes fl5|| . The 
same permutation sums (see eq. ( p. 13 )) and subsequent cancellation of tree- level amplitudes 



occurs in the collinear splitting of two-quark sub-leading partial amplitudes as well. 

The primitive decomposition of the splitting amplitudes is thus not important to the com- 
putation of the one- loop cross section in the collinear limit. This should not be too surprising 
since collinear factorization is a global property of the entire tree-level amplitude, not just of 
the color-ordered amplitudes, which play the role of tree- level primitive amplitudes. Nonethe- 
less, the primitive splitting rules are important for understanding the structure of primitive 
amplitudes and serve as useful checks on the primitive amplitudes themselves. Moreover, the 
primitive decomposition is important for an understanding of universal soft factorization. Even 
at tree-level, soft factorization is only a property of color ordered sub- amplitudes, not of the 
amplitude as a whole. Since soft factorization becomes apparent at the primitive level, it is 
convenient to express all infrared behavior at that level. 

We can derive the primitive splitting rules for the all-gluon and the two-quark amplitudes 
directly from the splitting rules of the leading partial amplitude. The reason for this is the 
simple primitive decomposition rule for the leading partial amplitudes. For instance, for n gluon 
amplitudes, 

A n:1 (l, 2, . . . , n) = 2, . . . , n) + ^-A% 2] (1, 2, . . . , n) + 2, . . . , n) , (4.42) 
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where A^\, J = 1, 1/2,0 are the primitive amplitudes. Corresponding to these primitive ampli- 
tudes are three one-loop primitive splitting amplitudes, Split '"^. The leading terms in N c in the 
splitting relation for A n .± are part of the primitive splitting function Split'- 1 -' , those proportional 
to nf/N c are part of Split' 1 ' 2 ', and those proportional to n s /N c are part of Split' '. That is, 



1112 



Split tre M„_ 1;1 + Split' 1 '^ + ^Split' 1 / 2 ^ 6 ! + ^Split'oU^ , (4.43) 



which breaks up into terms like eq. (|3,5|), 

M ^ SpUt^W + Split' 1 ^ , 



A m IE spiit^Mj^ + spiit' 1 / 2 '^ 



(4.44) 



A 



12 



nil 



SpUt* 166 ^!.! + Split^U^ , 



where the helicity sum has been suppressed. 

Similarly, when considering amplitudes with a pair of external quarks, there are four types of 



primitive amplitudes, An'^\ An'^\ A 1 ^ 1 ^ 2 ^ and An'^ in eq. ( 2.11 ), and, accordingly, four prim- 
itive splitting amplitudes. The leading terms in iV c in the splitting relation for A n -\ correspond 
to Split 1- , those suppressed by l/iV 2 correspond to Split R , etc. 



4.5.1 Primitive splitting rules for all-gluon amplitudes 



Using eq. fl4.44|) and the results of section O, the splitting rules for all-gluon primitive ampli 
tudes are 



Split" 9 ' [1] (1+2+) 



r" 9 ' [1 '(±,l T ,2*) 



ng,[l/2],, ,± 



+ 



rl 9, '^(A, 1 ± ,2 T ) 



Spht7 [1/2] (l + ,2+)= + 



r ^[V2] (±)lT)2T) 



(A,l ± ,2^) = 0, 



Sl2 t 
2 ' 



A* 



M \ [12] 2^z(l-z)(l-e5 R )cr 



si2 (12) 2 (l-2e)(2-2e)(3-2e) 



1 — z \ e ire 



z J sin(-7re) 



+ £ 2e 2m " 1 Li 2m _ 1 ( - 



m=l 



2z(l - z)(l - 5 R e) 



-s 12 (l-2e)(2-2e)(3-2e) 



1 — z \ e 7re 



z / sin(-7re) 



-«12 / 

fi 2 V [12] 2^z(l - z)c r 



m=l 



si2 (12) 2 (l-2e)(2-2e)(3-2e) ' 



2z(l - z) 



-si 2 i (l-2e)(2-2e)(3-2e) ' 



+ £ 2e 2m " 1 Li 2m _ 1 ( - 



(4.45) 



(4.46) 
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Split 



I* 



[12] 



2y/z{l - z)c r 



] (±,1 T ,2^)= + 



ft 



S12) (12) 2 (l-2e)(2-2e)(3-2e) ' 
2 \' 2z(l-z) 



-«ia / (l-2e)(2-2e)(3-2e) ' 



(4.47) 



r^' [Ol (A,l ± ,2^) = 0. 



4.5.2 Primitive splitting rules for two-quark amplitudes 

It is clear from eq. ( 2.1 1| ) that we can obtain splitting rules for at least some of the primitive 



amplitudes from the results of section O. However, there are additional primitive amplitudes 
for two-quark processes, those where the anti-quark and quark are not adjacent in the ordering 
like A§ (If) 3,2 g , 4, 5), that contribute to sub-leading, but not the leading, partial amplitudes. 
Nonetheless, it turns out that all of the rules can be derived from the splitting of the leading 
partial amplitudes. 

In g — > gg splitting, there are two possibilities: the collinear pair could lie between the quark 
and the anti-quark in the ordering, (i.e. on the gluonic part of the loop) as it must in the leading 
partial amplitudes, or, as it can in the additional primitive amplitudes, the collinear pair could 
lie between the anti-quark and the quark (on the fermionic part of the loop). Similarly, the 
leading partial amplitudes can have q — > qg or q — > gq splitting, but the additional primitive 
amplitudes can also have q — > gq or q — » qg splitting. Clearly, neither g — » qq nor g — > qq 
splitting contributes to the splitting of primitive amplitudes where the quark and anti-quark are 
not adjacent in the ordering. 

It would seem that we must take these extra cases into account. However, because the 



Ai s are related to the A n s by a reflection identity, eq. Q2.12 ), in which the ordering is 



reversed, an 'i?'-type primitive with the collinear pair lying between the quark and anti-quark 
is equivalent to an "L"-type primitive with the collinear pair lying between the anti-quark and 
the quark. The 'extra' splitting processes are thus already taken account in the splitting of the 
leading partial amplitudes. 

We thus need only to derive the splitting rules for the 'L'-type primitive amplitudes. There 
are eight cases: g — > gg splitting when the collinear pair lies between the quark and the anti- 
quark; g — > gg splitting when the collinear pair lies between the anti-quark and the quark; 
q — > qg splitting; q — > gq splitting; q — > gq splitting; q — > qg splitting; g — > qq splitting; and 
g — > qq splitting. 

For g — > <7ig2 splitting when the collinear pair is between the quark and anti-quark the 
splitting functions are identical to the pure glue case of eqs. ( |4.45 ), ( [4.46 ) and ( 4.47| ): 



Split^^T,A 1/2 '° ] (l Al , 2 Aa ) = Split™ 5 ' [^=1.1/2,0] (1 a 1;2 a 2) _ (448) 
For g — * g\g2 splitting when the collinear pair is between the anti-quark and quark, 

Split^i'^ 172 ' 01 (1 Al , 2 A2 ) = . (4.49) 
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-«12 



For g -»■ gig 2 splitting, 

2g,L,[l] 1± 2=F ^ _ 

29, Ml/2] Q ± T , _ 
25,L,[0] (x 1± 2T x _ 

For g ^ qiq 2 splitting, 

^•W,(A,l±,2 T ,= - 
r^] (A . 1± , 2T) = . 



+ 



13 



1 



1-z 



+ 



ire 



+ E 2e 2m - 1 Li 2m _ 1 f - 



sin(7re) 

v ' m=l 

5 - 2e - 3<5^ 



/i 



«12 y 
2 s 



6 e(l - 2e) 3(1 - 2e)(2 - 2e)(3 - 2e) 

2 1 2 

+ 



-S12 



3e(l-2e) 3(1 - 2e)(3 - 2e) 
1 1 2 



3e(l-2e) 3(1 - 2e)(3 - 2e) 



(4.50) 



A* 



-S12 



1 3 



1 l + <5 fl (l-2e) 



e 2 2e(l-2e) (l-2e)(2-2e) 



(4.51) 



In the following z refers to the momentum fraction of parton 1. Thus, in q — > (7152 splitting, z 
is the quark momentum fraction, while in q — > 51^2 splitting, z is the gluon momentum fraction. 
For q — > (7152 splitting, 



(±,1 T ,2 T ) 



29, L, [1] 

„2 9 ,Mi] (±)lT)2±) 

-A ±A 



A* 



1 



/X 



S12 / 

2 ~ 



' 5 

29, M1/2] / A , ± _ 
2?, MO] /x -,-A 2 ±n _ n 

For q^ giq 2 splitting, 



S12 



m=l V 

00 , 

l-^e m Li m f — 



(l-z)(l-e^) ) 
(l-2e)(2-2e) J ' 



m=l 



(4.52) 



29, Mi] 

S[q^gq] 

r 2g,Ml] 



29, L, [1/2], 
S [q^gq] 

2q, L, [0] 
r S[q^gq] 



(±,1 T ,2^) 



(±,l T ,2 d 



(A,1" A ,2 ± ) = 0, 



/i 



Sl2 / 
2 ' 



-«12 



]Te m Li m — 

m=l V 

00 ✓ 



z(l - egg) 1 
(l-2e)(2-2e)J ' 



m=l 



(4.53) 



(A,1" A ,2 ± ) = 0. 
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For q — ¥ gxq% splitting, 



2q,L, [1] 
S [q-igq] 

_2ff,i, [1] 



(±,1 T ,2^) 



S [q~>gq] 
2q,L, 



(±,l T ,2 d 
l (A,l" A ,2 d 



A ,2 ± ) 



' Slq-.gq]^ 1 

For q — ► qi52 splitting 



„2q,L, [1] 
„2q,L, [1] 



(±,1 T ,2^) 



(±,1 T ,2 ± ) 



S [q~*qg] 
r 2q, L, [0] 
S[q^qg] 



(A,1~ A ,2 ± ) = 



0, 
0. 



2q,L,[l/2] (x A 2 ±N _ 



Sl2 / 
2 ' 



"Sl2 



2 ' 



1-5: e - Li 

m=l 

1- 5> m Li m 



m=l 



1 - Z 



1 - 2 



*(1 ~ egg) 
(1 - 2e)(2 - 2e) 



oo 



-«12 / e 



m=l 



1 - Z 

—z 
— z 



:i-z)(l-e5 R ) 
(1 - 2e)(2 - 2e) 



(4.54) 



(4.55) 



As for the leading color partial amplitudes it is straightforward to rewrite the primitive 
splitting amplitudes in terms of formal polarization vectors and spinors. 



4.6 One-loop soft amplitudes 

The situation with soft limits is quite similar to the collinear case. Again the contributions 
can be divided into factorizing and non- factorizing pieces. The factorizing contributions to 
the soft limit of one-loop n point amplitudes come from those diagrams where an n — 1 point 
tree diagram emits a soft gluon via a loop. These are exactly the same diagrams that gave 
the factorizing contribution to the collinear splitting amplitudes, although we are now taking a 
different infrared limit. Taking the gluon momentum k 2 to zero in eqs. ( [4.2|) and ( [4.14 ), we find 
that the factorizing contributions to the one-loop soft amplitudes vanish. 

Just as the non-factorizing contributions to collinear splitting came from taking the limit of 
infrared singular loop diagrams, the non-factorizing contributions to the soft amplitudes come 
from taking the soft limit of those same loop diagrams. In table || we list the infrared divergences 
and the corresponding discontinuity functions that can appear in the soft limit. 

The complete one-loop soft amplitude as the momentum of leg 1 vanishes is given by 
5 1 - loop (n, 1, 2) = S iact (n, 1, 2) + Sl f 2 (s 12 , e) + s 2 f 2 (s nl ,e) + s 3 f 5 (s nl , s 12 , s n 2, e) , (4.56) 



where f 2 is defined in eq. ( 4.17 ) and /s is (see appendix C.2 ) 



/5(sni,si2,s n2 ,e) = lim . -fi A€ s nl s 12 l^ = --^ I ^ ( s n2 ) \ _ (457) 

fci^o2 e 2 \{-s n i)(-si2) I sm(7re) 
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Singularity 


Non-Factorizing Contribution 


c r ln(-s n i) 
e 


/20„i,e) 


crln(-si 2 ) 
e 


/2(si2,e) 


c r ln((-s„i)(-si 2 )/ (-s n2 )) 


/5(Snl)Sl2,Sri2,e) 


e 



Table 2: The potential 'non- factorizing' contributions to the soft amplitudes for ki — > 0. 



Since 5 fact (n, 1,2) =0, all that is needed to get the complete soft amplitude is to examine 
the infrared singularities of one-loop amplitudes and adjust the coefficients Sj so that we have 
the correct poles in e. 

Once again we can obtain the primitive soft rules for zero- and two-quark amplitudes from 
the soft rules of the leading partial amplitudes. Using the known divergence structure of the 
soft amplitudes ( 3.1C| ), we find the infrared divergence structure, 



-. non- fact 



--c r S 



tree 



1 ln(-s n i) ln(-si 2 ) ln(-s n2 ) 



--crS 



tree 



e z 

1 ln(- 

-^7 + — 



+ 



"Snl)(- 



e 

Sl2)/(- 



■Sn2 



+ O(e ) 



+ O(e ) 



(4.58) 



Using table g to match these divergences to their proper discontinuity function, it readily follows 
that 



-.l-loop / 



(n, V t , 2) = 5 tree (n, 1±, 2)/ 5 (s„i, s 12 , s n2 , e) 



%S tme (71,^,2) 



/^ 2 (-Sn2) 



ire 



(4.59) 



sin(-7re) 



v (-s„i)(-si 2 ) 

in agreement with our earlier result f2§. Although the helicity form of S tvee (n, 1 ± , 2) 

is usually the most convenient representation, within the context of conventional dimensional 
regularization it is convenient to express 5 tree (n, 1 ± , 2) in terms of the eikonal factor in eq. (ft .71). 

This result does not depend on 6r, so it is independent of the variety of dimensional regu- 
larization scheme used. This result is also independent of the particle and helicity types of the 
neighbors n and 2; a soft gluon gives rise to the same factorization term when both neighbors are 
gluons as it does when one or both neighbors are quarks. However, it does depend on whether 
the soft gluon is attached to a gluonic part of the loop or to a fermionic part of the loop. In 
eq. ( p. 11 ) all gluons in the leading in colors contribution are attached to the gluonic part of 
the loops, while in the sub-leading in color contributions they are attached to the fermionic 
(or scalar) parts of loops. Since there is no rif, n s or N c dependence in eq. ( 4.59| ) there are 
no sub-leading contributions. This property may also be expressed in terms of the primitive 
amplitudes, as we do below. 
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4.6.1 Soft factorization of all-gluon primitive amplitudes 



We see from eq. ( 4,59 ) that there are no nj or n s dependent contributions to soft factorization. 
We therefore have the simple result that 



?ng, [1] 



(11,1*2) 



S tree (n, 1*, 2) 



TTC 



H 2 {-Sn2) \ 

(-s n i)(-si2) / sin(7re) ' 



5™ 9 ' [1/2] (n,l ± ,2) = 0. 
S naM (n,l±,2) = 0. 



(4.60) 



4.6.2 Soft factorization of two-quark primitive amplitudes 

As in the case of primitive splitting relations, we must distinguish between situations where the 
soft gluon falls between the quark and the anti-quark in the ordering and those where it falls 
between the anti-quark and the quark. Again casting all factorization properties in terms of 
the 'L'-type primitive amplitudes, we find that when the soft gluon is between the quark and 
the anti-quark (i.e. the soft gluon is attached to a gluonic part of the loop), the one loop soft 
amplitudes are 

5 2 «' z 'W(n,l± 2)= -^< ? t-e (M±)2) /' ^~ s ^ V ™ 



e 2 \(-Sni)(-si2) J sin(7re) ' 

5 2 9 ,L,[i/ 2](n;1±)2) = 0; (4-61) 

S^-MO] l( n ,l±,2) = 0. 

When the soft gluon is between the anti-quark and the quark so that it attaches to a fermionic 
part of the loop, the one-loop soft amplitudes vanish, 

5 2<? ' L ' [1] (n,l ± ,2) = 0, 
5 2g.i,[l/2]( n> l± j 2) = 0, (4.62) 
s 2q,L,[0] ( n>1 ± i2 ) = 0. 

The case of four or more quarks is similar: if soft gluon in any primitive amplitude connects 
to a fermionic part of the loop then the soft factor vanishes. If the soft gluon is attached to a 
gluonic part of the loop then the soft factor is given by the universal function in eq. ( |4.61| ). 

4.6.3 Example: The soft limit of A 5;3 (%, 2 g ; 4, 5; 3). 

To illustrate the utility of the primitive decomposition of the one-loop soft and collinear splitting 
amplitudes, we compute the soft limit of one of the sub-leading partial amplitudes for two-quark 
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three-gluon scattering. Using eq. ( |2,15 ), and taking the limit that gluon 3 becomes soft, 
A5. 3 (l ff ,2,;4,5;3) 



A%(l q ,4,5,2 q )+Ai(lq,5,4,2 q ) 



-.tree , 



1,2) 



+ A%(l q , 2 q , 4, 5) 5 tree (2, 4) + S tree (4, 5) + S tree (5, 1 



+ ^(1 ? ,2 9) 5,4) 
+ A%(l q ,4,2 q ,5) 



5 tree (2, 5) + 5 tree (5, 4) + 5 tree (4, 1 



S tree (2,5)+S trcc (5,l) 
5 tree (4,l) +5 trec (2,4) 



(4.63) 



< ce (l 5 , 2 q , 4, 5) [5 1_loop (l, 2) + 5 1 " loop (2, 5) + S 1 " 100 ^, 4) + 5 1 " loop (4, 1) 

(1, 2) + 5 1 " loop (2, 4) + 5 1 " loop (4, 5) + S 1 ' 10 ^^, I) 



A 



| ree (l g -,2 ? ,5,4) 



-.1-loop 



where we abbreviate S tiee (i,j) = S tiee (i,3,j) and 5 1 " loop («',i) = S 1 ' 1 ^^, 3, j) 
properties for taking the limit are 5 tree (i, 3, j) = — S tTce (j, 3, i) and 



Some useful 



Af e (l g -,4,2 9 ,5) =Ar e (l q ,5,2 q ,A) = -A^ e (1 ?> 2 q , 4, 5) - < ee (l g , 2„ 5, 4) , 



(4.64) 



As mentioned in section |2| the tree-level primitive amplitudes correspond to the partial ampli- 
tudes when all particles including the fermions are taken to be in the adjoint representation. 
These properties of the tree-level primitive amplitudes are simply the properties Q of any such 
partial amplitudes. 

Comparing this result to eqs. ( [2.14 ) and (2.15), it is clear that there is no simple soft factor- 
ization relation for sub-leading color partial amplitudes. In particular, there is no direct analog 
to eq. ([O]) which would express the factorization in terms of lower point partial amplitudes. 
The reason for this is that soft factorization depends upon the soft gluon's neighbors in the color 
flow. Partial amplitudes, especially sub-leading partial amplitudes, can have a rather compli- 
cated color flow [24] with the result that soft factorization becomes deeply entangled. Primitive 
amplitudes have a fixed ordering of external legs so there can be no twist in the color flow to 
complicate soft factorization. By piecing together these primitive components, the soft limits of 
sub-leading partial amplitudes can be systematically understood. 



5 Checks on results 

5.1 Checks using supersymmetry. 

Although QCD is not a supersymmetric theory we may use supersymmetry as a way of verifying 
our results. As described in the second appendix of ref. [^], we may convert the QCD splitting 
functions to those for N = 1 super- Yang-Mills theory, which contains a single gluon and gluino, 
by adjusting the color factors to be — ► —1, nt/N c — ► 1 and n s /N c — > in the splitting 
amplitudes (eqs. ( 4.32 ) - ( 4.37] )). (As discussed in appendix B of ref. (2f|, one can also check 



the case of an iV = 1 chiral multiplet, but then one must account for Yukawa interactions.) 
To preserve the supersymmetry we take the parameter controlling the variant of dimensional 
regularization to be Sr = 0. 



27 



With these substitutions, the only non- vanishing splitting amplitudes are those where the 
tree-level amplitudes do not vanish. All the remaining non-vanishing splitting amplitudes are 
given by 

Split_ A (l Al ,2 A2 ) =c r Split tr 5 l c (l Al ,2 A2 )r 5 (-A,l Al ,2 A2 ), (5.1) 

where 

r s (-X, l Al , 2 Aa ) = f 2 (s 12 , e) + / 3 (a 12 , e, z) + f 3 (s 12 , e,l-z), (5.2) 

is independent of the specific helicities or particle labels. (The explicit value of the sum of these 
functions is given in eq. ( C.iop .) This independence is a consequence of supersymmetry p3] as 



has been previously noted through O(e ) [26]. The fact that our expressions (for 5r = 0) satisfy 
this constraint to all orders in e provides non-trivial support that the results are correct. 

For the special case of N = 4 super- Yang-Mills explicit results for the one-loop four- five- and 



six-point amplitudes are known to all orders in e [34]. Using these results, we may extract the 
splitting amplitudes simply by taking the collinear limit of the five-point amplitude and com- 
paring it to the four-point amplitudes. The result so obtained agrees with the supersymmetric 



results in eqs. (5.1) and fl5.2|) , providing another independent check. 

We have also performed the same checks in the case where an external gluon momentum 
becomes soft. In this case the supersymmetry is manifest since the soft amplitude does not 
depend on the number of flavors or colors or on whether the neighboring particles are gluons or 
fermions. 



5.2 Checks from Higgs amplitudes 

We have also checked the one- loop g — > gg and g — > qq splitting amplitudes, as well as the 
one-loop soft amplitude, by directly computing to all orders in e the one-loop gggH and gqqH 
amplitudes in a theory with an effective ggH coupling [35| 



£eff = ~\ 



a q H ( 11a, 
1 -— 1 + 



37r v \ An 



Tr G^CT, (5.3) 



given by the infinite mass limit of a heavy quark triangle. The gggH and gqqH amplitudes in 
this theory are convenient for extracting the splitting and soft amplitudes since they involve only 
four-point kinematics with a single massive leg; the massive Higgs ensures that both amplitudes 
have well-defined soft and collinear limits. Since there is no tree-level qqH coupling in this 
effective theory, we cannot use these processes, however, to check q — > qg splitting. 

The details of the calculation through (D(e ) may be found in ref. |^6|. For this comparison 
we have extended the calculation to all orders in e. We find that the amplitudes can be written 
as simple analytic expressions with the scalar box integrals expressed in terms of hyper geometric 
functions, as given in appendix [C]. 

Since there are only three external partons, there are no sub-leading partial amplitudes for 
these processes. Primitive decompositions are trivial to perform. The gggH amplitudes mirror 
the all-gluon amplitudes of pure QCD with the primitive amplitudes given simply by the parts of 
the leading partial amplitude corresponding to the different spins of the loop particle. The gqqH 
amplitudes mirror the two-quark amplitudes of pure QCD, with the rules for those orderings 
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with the gluon between the quark and anti-quark (lg,2 g ,3) differing from the rules for those 
with the gluon between the anti-quark and quark (1^, 3, 2 q ). 

With the gggH amplitudes, we have checked the one-loop g — > gg splitting amplitudes 
(eqs. ( 4.32; ) and ( [4.33j )) as well as the one-loop soft amplitudes (eq. ( |4.59| )). With the gqqH 
amplitudes, we can check the one-loop g — > qq splitting amplitudes (eqs. ( |4.35 ), ( 4.50 ) and 
( |4.51 )). In all cases, we find complete agreement to all orders in e. These checks emphasize the 
universal, i.e. process independent, nature of infrared factorization. 



6 Conclusions 

In this paper we have presented the functions describing the behavior of one-loop amplitudes 
in the infrared divergent regions of phase space. These functions should be useful for NNLO 
n-jet computations. When evaluating the phase space integrals associated with the (n + 1)- 
point one-loop amplitudes with one unresolved parton, infrared divergences of order e -2 in the 
dimensional regularization parameter are encountered. This suggests that one needs to know 
the one-loop amplitudes through 0(e 2 ). Given the complex analytic structure of multi-parton 
amplitudes || ||, [2^, 23, 24, 25], it would be rather non-trivial to obtain the complete higher 



order in e contributions. A better approach would be to use the splitting functions presented in 
this paper to all orders in e to obtain the higher order in e contributions only in the relevant soft 
and collinear regions of phase space. Then the only calculations needed would be the (n + 1)- 
point one-loop amplitudes through O(e ) and the n-point amplitudes through C(e 2 ), which is 
much simpler to obtain. We also discussed the decomposition of QCD amplitudes into primitive 
amplitudes as a way to disentangle the color factors from the infrared divergent regions of phase 
space; for the case of external fermions these primitive amplitudes are a finer subdivision of the 
amplitudes than the conventional color decomposed ones. 

In order to extend the analysis in this paper to the case of two-loop soft and collinear 
splitting amplitudes, one would need a detailed understanding of the infrared structure of two- 
loop integrals as well as of the discontinuity functions which may arise in the soft or collinear 
limits. The factorizing contributions are, however, straightforward to compute since they involve 
only triangle contributions. 

The development of formalisms that would enable one to compute the NNLO contributions 
to multi-jet amplitudes remains an important open problem in perturbative QCD. There has, 
however, been some recent progress towards this. For example, the singularity structure when 
two partons are unresolved has been examined in ref. |2C|| . An interesting formalism for obtaining 
collinear splitting amplitudes from unitarity |2[ 26, ^] has also recently been developed 



Furthermore, two-loop four-gluon amplitudes in the special case of maximal supersymmetry 



have been evaluated in terms of scalar integral functions [12] using the new techniques described 
in ref. 0]. Although a number of substantial difficulties remain, these developments provide for 
some optimism. 
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A Tree-level splitting amplitudes 

We collect here the helicity form of the tree- level splitting amplitudes GJ, using the sign con- 



ventions of ref. [24|. For the case of a gluon splitting into two gluons, with helicities labeled as 



if all particles are outgoing, the splitting amplitudes are 

,g-»SlS2,tree ri + +x _ n n luP^Si^, tree/- - - ~1 



Split™ 2 ' troc (l + ,2+) =0, s P ht™ 2 ' tree (i-,2-) 

S P ht™ trcc (l-,2+) = ^(1-%(12) ' Spht™ trcC (l + ,2-) 
Sp ttt™ 2 ' tree (l-,2-) = 0, Split^ 9192 ' trcc (l + ,2+) 



Spllt - (1 ' 2 >- y/z(l-z)[12] ' Spht™ 2 ' tree (l-,2^ 
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(A.l) 

The spinor inner products are (ij) = and [ij] = where \i^) are massless 

Weyl spinors of momentum fcj, labeled with the sign of the helicity. They are antisymmetric, 
with norm | (i j) \ = \ [ij] \ = y / s7J, and obey the rule (ij) [j i] = Sy = 2ki ■ kj. 

For g -> gig 2 , 

z l / 2 (l - z) 3 / 2 - z 3 l 2 (\ - z) 1 / 2 

QnlH^'?i<32,tree/- | + -\ _ f I 1 D„m-9-»9l92, tree /-, - _ f V 1 f ) 
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- * z 1 / 2 (l - z) 3 / 2 - * z 3 l 2 (\-z\ 1 l 2 

q„iu9-*9l92,tree/-,- 9 +\ _ f I 1 f J c r .i: + S-»9i92,tree/ 1 + 9 _\ _ f I 1 

Split. ( i_,2 ? )--^=^, Split. ( i_, 2? 

(A.2) 
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Finally, for q -> q^ 2 and q -> gi<j 2 , 
Splitr 9192,tree (l+,2+) = 

Split™ 32 ' tree (l-,2-) = 
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(A.3) 

Because we are working with color ordered amplitudes, we must distinguish between, say, 
q — > q\gi and q — ► 51(72- This is accomplished by exchanging 1 and 2 (and therefore 2 and 1 — z). 
For example. 



- ♦ ^3/2/-, _ „\l/2 3/2(1 \l/2 

p + ( 9 ' ? j " yia^y (2 1) " v^a^y u 2) 



(A.4) 



B Previous results through O(e ) for one-loop splitting ampli- 
tudes. 

For convenience we collect the splitting and soft amplitudes through 0(e). These were obtained 
from refs. (2^, |24}| . These are useful since they contain the infrared singularities which are used 
to fix the non-factorizing contributions. They also provides a check on our results. 

The loop splitting amplitudes have a structure similar to the tree splitting amplitudes, so it 
is useful to express them in terms of a proportionality constant rs defined by 

Split ^ oop (l Al ,2 A2 ) = c r x Splits (l Al ,2 Aa ) x r s (-A, 1 A \ 2 Aa ) , (B.l) 

for general partons 1 and 2. The only exception to eq. (B.l) is for g~ — > g + g + (and its 



parity conjugate g + — ► g~g~), where Split tree vanishes but Split 1 " 10013 does not. In general 
rg(— A, l Al , 2 A2 ) depends on the parton helicities. As a notational point we distinguish between 
particles circulating in the loop by a spin index, i.e. Split [J] (1, 2). For the cases J = 0, 1/2, 1 
respectively correspond to a scalar, a fermion, and a gluon in the loop. 

The Splitt^ (1 + , 2 + ) obey the supersymmetry relation Split^ = — Split^ 1 / 2 ' = Split ^, where 
SplitW(l+,2+) = -^z(l-z) . (B.2) 

We present the remaining g — ► gg loop splitting amplitudes in terms of rs- 

r ™i-loo P(±) 1T)2T) = _l_ (l* V + 2 In. ln(l -z)~ 

e z \z(l - z)(-si2) J 6 

1 ,„ « /„ ,. nt n s 



rr 9lfl2,1 - l00P (A, 1*7?) = -~ ( ^— C \ J + 2\nz ln(l -*)"£■ 

e 2 \z(l - z)(-si 2 ) I 6 
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The rs functions for the case 
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C Evaluation of box integrals in collinear and soft limits 

In this appendix we evaluate the integral functions appearing in the one-loop soft and collinear 
splitting amplitudes to all order in e. In the collinear or soft limits, the hypergeometric functions 
in terms of which the integrals may be expressed reduce to sums of poly logarithms. These results 
allow us to present explicit formulas for the non-factorizing contributions, valid to all orders in 
dimensional regularization. 




Figure 5: The particular box integral under consideration in this appendix. Leg 3 is composed of a 
sum over massless momenta so that K\ ^ 0. 
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Consider the single-mass box integral Zj.g^s, t, iff ), depicted in fig. ||, where we take the 
d leg to be massive, iff 7^ and the Mandelstam variables appearing as arguments are, 

s = {k 1 + k 2 f t={k 1 + k 4 f. (C.l) 



Limits of this integral appear in the collinear functions fa in eq. ( [4.17| ) and in the soft function 
in eq. ( |C.13[ ). All other box integrals appearing in the soft and collinear splitting amplitudes 
are given by relabelings of this integral. 

It can be expressed as a sum of hypergeometric functions JI^ ] 
2c r n 2e f ft- K f 



-ifi 2e stll% = i \ -— ^ I 2 FA -e, -e; 1 - e; 1 + 



e 2 \\ st J z \ ' ' ' ' t-K$ 
+ ( S - L ^) 6 2F 1 (-e, -e; 1 - e; 1 + — L,) (C.2) 



( S -iff)(t-iff)^ f , F , *« 

2-^1 — e, — e, 1 - e, 1 



C.l The collinear limit of the single-mass box integral 

In the collinear limit k\ = zP and k 2 = (1 — z)P, we find that s — > 0, t — > ziff, and eq. (|C. 
reduces to 

- iip8tZl£(8, t, iff)* ^ ^ 2 F X (-6, -e; 1 - e; 1 - z) . (C.3) 
We wish to re-express this as a power series in e. Using the hypergeometric identity 

2 F 1 (-e, -e; 1 - e; 1 - z) = z' 2 F X (l, -e; 1 - e; "^^) , (C4) 
and the expansion of the hypergeometric function as a power series in z, 



- (i_ e )( 2 _ e )... („_£)(„.) 

00 / 00 n \ 00 

!-E^ m Ei = i-E^ 



(C.5) 



m=l \n=l / m=l 



we can write the single external mass box integral (eq. ( |C,3| )) as 



(C.6) 



Note that this contains the infrared singularity, 



2cr T . / 2cr , , ,_, 

— L11 f- v - M = — ln(z) . (C.7) 
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As discussed in section Q the known coefficient of this singularity may be used to determine 



the overall coefficient of the function fs(si2, e, z) given in eq. ( 4.17 ). Similarly, the coefficient of 
fz(s\2, e, 1 — z) is fixed from the coefficient of 2cr ln(l — z)/e. 

When the splitting amplitude contains the combination fy(s, e, z) + fz{s, e, 1 — z), as in the 
pure gluon case, we may further simplify the expression. Using hypergeometric identities (see 
eq. (9.132) of ref. §1), we find that 



2*1 (l, 
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= 1-2*1 (1, e-l + e-z) + (-js)-T(I - e)r(l + e) 
so that we can rewrite —iii 2tL stT\™(s : t) as 
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Thus, using eq. d4T7| ), the total non-factorizing contribution for the case of g —* gg splitting, 
given in eq. (|4.26| ) is 

rs = f2(s, e) + f 3 (s, e, z) + f 3 (s, e, 1 - z) 
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z ) sin(7re) 

C.2 The soft limit of the single- mass box integral 



In the soft limit k\ — * 0, we have s — > 0, t — > 0, and the single mass box integral eq. (|C.2; ), 
reduces to 



2cr ( vH-Kjy 
e 2 \(s)(-t) 



J F 1 (-e,-e;l-e;l) . 



Using the identity 



l F 1 (-e, -e; 1 - e; 1) 



7Tf 
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we can write the soft discontinuity function f5(s n \, s\2, s n2 , e) to all orders in e as 
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